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The path-integral quantization of thermal scalar, vector and spinor fields is performed newly 
in the coherent-state representation. In doing this, we choose the thermal electrodynamics and ip 4 
theory as examples. By this quantization, correct expressions of the partition functions and the 
generating functionals for the quantum thermal electrodynamics and (p 4 theory are obtained in the 
coherent-state representation. These expressions allow us to perform analytical calculations of the 
partition functions and generating functionals and therefore are useful in practical applications. 
\f} • Especially, the perturbative expansions of the generating functionals are derived specifically by 

f*"^ ' virtue of the stationary-phase method. The generating functionals formulated in the position space 

, are re-derived from the ones given in the coherent-state representation. 
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I. INTRODUCTION 



In the quantum statistics, it has been shown that the partition functions and thermal Green's functions for many- 
body systems may conveniently be calculated in the coherent-state representation [1-6]. This is because the partition 
• function either for a boson system or for a fermion system can be given a path (or say, functional) integral expression 
\ in the coherent-state representation and, furthermore, one can write out a generating functional of Green's functions 
■ for the system in the coherent-state representation [1-6]. However, the path integral expression derived in the previous 
literature can only be viewed as a formal symbolism because in practical calculations, one has to return to its original 
discretized form which leads to the path integral expression. If one tries to perform an analytical calculation of the 
path integrals by employing the general methods and formulas of computing functional integrals, one would get a 
wrong result. This implies that the previous path integral expressions for the partition function and the generating 
functional of Green's functions were not given correctly. The incorrectness is due to that in the previous path integral 
expressions, the integral representing the trace is not separated out and the time-dependence of the integrand in the 
remaining part of the path integral is given incorrectly. The partition function and the generating functional of Green's 
functions were rederived in the coherent-state representation and given correct functional-integral expressions in the 
author's previous paper [7]. These expressions are consistent with the corresponding coherent-state representations of 
the transition amplitude and the generating functional in the zero-temperature quantum theory [9-11]. Particularly, 
when the functional integral is of Gaussian type, the partition function and the generating functional can exactly 
be calculated by means of the stationary-phase method without any uncertainty [9-11]. For the case of interacting 
systems, the partition function and finite-temperature Green's functions can be conveniently calculated from the 
generating functional by the perturbation method. 

The aim of this paper is to formulate the quantization of thermal fields in the coherent-state representation. To 
be definite, we will choose the thermal quantum electrodynamics (QED) and the thermal ip 4 theory [5,6] as examples 
to describe the quantization of scalar fields, fermion fields and gauge fields. In this quantization, we will first derive 
correct path-integral expressions of the partition functions and the generating functionals of Green's functions for 
these fields. Then, we focus our attention to the perturbation method of calculating the partition functions and the 
generating functionals. In the zero-order approximation, the partition functions and the generating functionals will 
be exactly calculated by means of the stationary-phase method, giving results as the same as those given previously 
from the theories formulated in the position space. 

The remainder of this paper is arranged as follows. In Sect. 2, we quote the main results given in our previous paper 
for quantum statistical mechanics. These results may straightforwardly be extended to the quantum field theory In 
Sect. 3, we describe the coherent-state representation of Hamiltonians and actions for the thermal QED and (p 4 
theory which are needed for quantizing the theories in the coherent-state representation. In Sect. 4, the quantizations 
of the thermal QED and (p 4 theory are respectively performed in the coherent-state representation. In doing this, 
we write out explicitly the generating functionals of thermal Green's functions for the theories mentioned above. In 
particular, we pay our main attention to deriving the perturbative expansions of the generating functionals in the 
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coherent-state representation. In Sect. 5, the perturbative expansions of the generating functionals given in Sect. 4 
will be transformed to the corresponding ones represented in the position space. In the last section, some concluding 
remarks will be made. 



II. PATH INTEGRAL FORMULATION OF QUANTUM STATISTICS IN THE COHERENT-STATE 

REPRESENTATION 

First, we start from the partition function for a grand canonical ensemble which usually is written in the form [4-6] 

Z = Tre-P* (1) 

where /? = with k and T being the Boltzmann constant and the temperature and 

K = H-nN (2) 

here \i is the chemical potential, H and N are the Hamiltonian and particle-number operators respectively. In the 
coherent-state representation, the trace in Eq. (1) will be represented by an integral over the coherent states. To 
determine the concrete form of the integral, it is convenient to start from an one-dimensional system. Its partition 
function given in the particle-number representation is 

oo 

Z = J2(n\e-^\n). (3) 
Then, we use the completeness relation of the coherent states [3-10] 



/ 



D{a*a) | a >< a* |= 1 (4) 



where | a > denotes a normalized coherent state, i.e., the eigenstate of the annihilation operator a with a complex 
eigenvalue a [9-12] 

a \a) = a \a) (5) 

whose Hermitian conjugate is 

(a*\a + = a* (a*\ (6) 
and D(a*a) symbolizes the integration measure defined by 

D(a*a) — \ \da*da, for bosons; 

^ ' da*da, for fcrmions. ^ ' 

In the above, we have used the eigenvalues a and o* to designate the eigenstates \a) and (a*\, respectively. Inserting 
Eq. (4) into Eq. (3), we have 

oo 

Z = ^2 D(a*a)D(a'*a') (n | a') (a'*| e~ pk \a) (a* | n) (8) 



where 



are the energy eigenf unctions given in the coherent-state representation (Note: for fermions, n = 0, 1) [4-10]. The 

both eigenfunctions commute with the matrix element ^a'* | | aj because the operator K(a + ,a) generally is 

a polynomial of the operator a + a for fermion systems. In view of the expressions in Eq. (9) and the commutation 
relation [3-10] 
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a* a' = ±a'a* (10) 

where the signs "+" and "— " are attributed to bosons and fermions respectively, it is easy to see 

(n | a'} (a* \ n) = (±a* | n) (n \ a'} . (11) 

Substituting Eq. (11) in Eq. (8) and applying the completeness relations for the particle-number and coherent states, 
one may find 

Z = J D(a*a) (±a* | e~ fik \a) (12) 

where the plus and minus signs in front of a* belong to bosons and fermions respectively. 

To evaluate the matrix element in Eq. (12), we may, as usual, divide the "time" interval [0,(3] into n equal and 
infinitesimal parts, (3 — ne. and then insert a completeness relation shown in Eq. (4) at each dividing point. In this 
way, Eq. (12) may be represented as [1, 3-6] 

n-l 

Z = J D(a*a) n D(a* ai ) (±a* \ e^ K \a n ^) 

i=i U-jJ 

x(al_ 1 \e- ek \a n - 2 ) (at\ e~ ek \a) . 

Since e is infinitesimal, we may write 

e - £i ^ a+ > a ) = l-eK(a*;a) (14) 

where K(a + ,a) is assumed to be normally ordered. Noticing this fact, when applying the equations (5) and (6) and 
the inner product of two coherent states [1, 3-12] 

(a* | Oi_i) = e -3«:°i-i°7-i««-i+»:<*-i (15) 

which suits to the both of bosons and fermions, one can get from Eq. (13) that 

n— 1 n 

Z = jD(a*a)e- a ' a J ft D(a*a,i) exp{-£ X) K ( a *, a i-i) 

(16) 

+ - a * a i} 

i=l i=l 

where we have set 

±a* = a* , a = a . (17) 

It is noted that the factor e~ a a in the first integrand comes from the matrix elements (±a* \ a„_i) and {al \ a) and the 
last sum in the above exponent is obtained by summing up the common terms and — ^a*_^ai-i appearing in 

the exponents of the matrix element (a* \ a^-i) and its adjacent ones (a* +1 \ di) and (a*_ 1 \ cn-2)- As will be seen in 
Eq. (21), such a summation is essential to give a correct time-dependence of the functional integrand in the partition 
function. The last two sums in Eq. (16) can be rewritten in the form 

n n— 1 

E a* a t -i - E a * a i 

1=1 n-l T 1 . (18) 
i=l 

Upon substituting Eq. (18) in Eq. (16) and taking the limit e — > 0, we obtain the path-integral expression of the 
partition functions as follows: 

Z = j D(a*a)e- a ' a j T){a*a)e^ a *^ (19) 

where 
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Q ^da* (T)da(r), for bosons; 
®( a a ) — { Y\da*(T)da(T), for fermions 

r 

and 

I (a*, a) = \a*(fi)a{P) + ±a*(0)a(0) - jj dr[\a* (r)d(r) 

-|d*(r)a(r)+A-( a *(T),o(r))] (21) 
= a*(/3)a(0) - dr[a* {T)a{T) + K (a* (t) , a(r)) 

where the last equality is obtained from the first one by a partial integration. In accordance with the definition given 
in Eq. (17), we see, the path-integral is subject to the following boundary conditions 

a*(/3) = ±o*,a(0) = a (22) 

where the signs "+" and "— " are written respectively for bosons and fermions. Here it is noted that Eq. (22) does 
not implies a(/3) — ±o and a*(0) = a*. Actually, we have no such boundary conditions. 

For the systems with many degrees of freedom, the functional-integral representation of the partition functions may 
directly be written out from the results given in Eqs. (19) -(22) as long as the eigenvalues a and a* are understood 
as column matrices a = (ai, 02, ■ ■ ■ , fflfe, ■ ■ •) and a* — (a*, a%, ■ ■ ■ , a k , • ■ •). Written explicitly, we have 

Z = J D(a* a)e- a >" J S(a*a)e /(a *' a) (23) 



where 



V[ ^da* k da,k , for bosons; 

D(a a) { Y\da* k da k , for fermions, 

k 

Y[ -da* k {r)dak{T) , for bosons; 

D(a*a) = { tt (25) 

v ' L \\da* k {T)dak{T) , for fermions v ' 

kr 



and 



I(a*,a) = a* k (f3)a k (f3)- / rfr[4(r)d fe (r) + K(a* k (r), a k (r))}. 
Jo 



(26) 



The boundary conditions in Eq. (22) now become 

aUP)=±at,a k (0) = a k . (27) 

In Eqs. (23) and (26), the repeated indices imply the summations over k. If the k stands for a continuous index as 
in the case of quantum field theory, the summations will be replaced by integrations over k. 

It should be pointed out that in the previous derivation of the coherent-state representation of the partition func- 
tions, the authors did not use the expressions given in Eqs. (16) and (18). Instead, the matrix element in Eq. (15) 
was directly chosen to be the starting point and recast in the form [1-6] 

(a* I = exp{-|[a*(^^-) - (l^Li )a ._ l]} . (28) 
2 e e 

Substituting the above expression into Eq. (13) and taking the limit e — ► 0, it follows [1-6] 



Z = J £(a*a)cxp{- J dr[^a* (t)ci(t) - ^d*(r)a(r) + K(a*(r), o(t))]}. 



(29) 



Clearly, in the above derivation, the common terms appearing in the exponents of adjacent matrix elements were not 
combined together. As a result, the time-dependence of the integrand in Eq. (29) could not be given correctly. In 
comparison with the previous result shown in Eq. (29), the expression written in Eqs. (19)-(21) has two functional 
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integrals. The first integral which represents the trace in Eq. (1) is absent in Eq. (29). The second integral is defined 
as the same as the integral in Eq. (29); but the integrand are different from each other. In Eq. (19), there occur two 
additional factors in the integrand : one is e~ a a which comes from the initial and final states in Eq. (13), another 
is e^ a (P) a (P)+ a (°) a (°)l in which a* {(3) and a(0) are related to the boundary conditions shown in Eq. (22). These 
additional factors are also absent in Eq. (29). As will be seen soon later, the occurrence of these factors in the 
functional-integral expression is essential to give correct calculated results. 

To demonstrate the correctness of the expression given in Eqs. (23)-(27), let us compute the partition function for 
the system whose Hamiltonian is of harmonic oscillator-type as we meet in the cases of ideal gases and free fields. In 
this case 

K(a*a) — u k a* k a k (30) 
where u) k = £k — M with e k being the particle energy and therefore Eq. (26) becomes 

I(a*,a) = aU(3)a k (p)- dT[a* k (r)d k (T) + ui k a* k { T )a k {T)]. (31) 
Jo 

By the stationary-phase method which is established based on the property of the Gaussian integral that the integral 
is equal to the extremum of its integrand [8-11], we may write 

S>(a*a)e J < * >0 > = e Io{a ' » (32) 

where Iq(cl* , a) is obtained from /(a*, a) by replacing the variables <z£(t) and a k (r) in I (a* , a) with those values which 
are determined from the stationary condition SI (a* , a) = 0. From this condition and the boundary conditions in Eq. 
(27) which implies Sa* k (fl) — and 5a k (0) = 0, it is easy to derive the following equations of motion [8-11] 

dfc(r) + Wfea fe (r) = 0, a* k {r) - Wfea^r) = 0. (33) 

Their solutions satisfying the boundary condition are 

a k (r) = a k e— T , a* fe (r) - ±a* k e^ T ^K (34) 

On substituting the above solutions into Eq. (31), we obtain 

I (a*,a) ^ ±a* k a k er UJkfi . (35) 

With the functional integral given in Eqs. (32) and (35), the partition functions in Eq. (23) become 

z _ , / D(a*a)e- a >*( 1 - e ~' 3 ^ , for bosons; . . 

° ~ 1 j D(a*a)e-< a " ( - 1 + e ' fiu ' k ^ , for fermions. ( ' 

For the boson case, the above integral can directly be calculated by employing the integration formula [1]: 



/ 



D(a*a)e- a ' (Xa - b \f(a) = jfiX^b). (37) 



The result is well-known, as shown in the following [3-6] 

z ° = Ilr^- (38) 

For the fermion case, by using the property of Grassmann algebra and the integration formulas [3-10] 

J da = j da* =0 , J da*a* = j daa = 1, (39) 

one may easily compute the integral in Eq. (36) and gets the familiar result [3-6] 

Z = l[(l + e~ f3 ^). (40) 
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It is noted that if the stationary-phase method is applied to the functional integral in Eq. (29), one could not get the 
results as written in Eqs. (38) and (40), showing the incorrectness of the previous functional-integral representation 
for the partition functions. 

Now let us turn to discuss the general case where the Hamiltonian can be split into a free part and an interaction 
part. Correspondingly, we can write 

K(a*,a) = K (a*,a)+H I (a*,a) (41) 

where Ko(a*,a) is the same as given in Eq. (30). In this case, to evaluate the partition function, it is convenient to 
define a generating functional through introducing external sources j k (r) an d ife( T ) sucn that [4,5] 

Z\j* j] = J D(a*a)e~ a " kak J D(a* a)e ak( -^ ak( - 13 ^^ dT i a k( T )'k(' r )+ K ( a * a )-fk( T ) a k(T)-a' k (T)j k (T)] 

= e -^ d ^*ifr>±^ Zolj *j ] (42) 

where the sings "+" and "— " in front of S j^( T ) TC ^ eT to bosons and fcrmions respectively and Z [j*,j] is defined by 

Z \j*,j] - J D(a*a)e- a ' kak J 3>(a*a)e J ( *> 0i ^'> (43) 

in which 

I(a*,a;j*,j) = al{/3)a k (P)-fodT[al(T)a k {T) (44) 
+cj k a* k (r)a fe (r) - j* k (r)a fe (r) - a* k (r) j k (r)] . 

Obviously, the integral in Eq. (43) is of Gaussian-type. Therefore, it can be calculated by means of the stationary- 
phase method as will be shown in detail in Sect. 4. 

The exact partition functions can be obtained from the generating functional in Eq. (42) by setting the external 
sources to be zero 

Z = Z\j*,j]\ j ,= j = . (45) 

In particular, the generating functional is much useful to compute the finite-temperature Green's functions. For 
simplicity, we take the two-point Green's function as an example to show this point. In many-body theory, the 
Green's function usually is defined in the operator formalism by [4,13] 

G k i(r u r 2 ) = iTr{e-^T[a k (ri)a+(r 2 )}} = Tr{e^' R )T[a k { Tl )a+ (r 2 )]} (46) 

where < n, T2 < /3, O = — ^ In Z is the grand canonical potential, T denotes the "time" ordering operator, a k (ri) 

and af(T2) represent the annihilation and creation operators respectively. According to the procedure described in 
Eqs. (12)-(22). it is clear to see that when taking t\ and Ti at two dividing points and applying the equations (5) and 
(6), the Green's function may be expressed as a functional integral in the coherent-state representation as follows: 

G k i(T U T 2 ) = \l D{a*a)e- a > k j ^(a* a)^^ (t^' ^ . (47) 

With the aid of the generating functional defined in Eq. (42), the above Green's function may be represented as 

Gfez(ri,r 2 ) = j«=j=o (48) 

Z 0] k {Tl)bjl{T 2 ) 

where the sings "+" and "— " belong to bosons and fermions respectively. 



III. COHERENT-STATE REPRESENTATION OF THERMAL FIELD HAMILTONIANS AND ACTIONS 

To write out explicitly a path-integral expression of a thermal field in the coherent-state representation, we first need 
to formulate the field in the coherent-state representation, namely, to give exact expressions of the field Hamiltonian 
and action in the coherent-state representation. For this purpose, we only need to work with the classical fields by 
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using some skilful treatments. In this section, we limit oursclf to describe the coherent-state representations of the 
thermal QED and <p 4 theory . 

Let us first start from the effective Lagrangian density of QED which appears in the path-integral for the zero- 
temperature QED [9,10] 

C = ${i-f(dp - ieA^) - M}yj - -F^F^ - ^-(0M„) 2 - d^Cd^C (49) 

4 la 

where and -tp represent the fermion field, is the vector potential of photon field, C and C designate the ghost 
field 

= d^Ay - d v A^ (50) 

and M denotes the fermion mass. It is noted here that although there is no coupling between the ghost field and the 
photon field, as will be shown later, it is necessary to keep the ghost term in Eq. (49). For the sake of simplicity, 
we will work in the Feynman gauge. In this gauge, the Lagrangian, which is obtained from the above Lagrangian by 
applying the Lorentz condition d^A^ = 0, is of the form 

C = ${i-f(dp - ieA„) - M}ip - ^A^A" - d^Cd^C (51) 

The above Lagrangian is written in the Minkowski metric where the 7— matrix is defined as 70 = (3 and 7 = [3d [10]. 
In the following, it is convenient to represent the Lagrangian in the Euclidean metric with the imaginary time t = it 
where t is the real time. 

Since the path-integral in Eq. (42) is established in the first order (or Hamiltonian) formalism, to perform the 
path-integral quantization of the thermal QED in the coherent-state representation, we need to recast the above 
Lagrangian in the first order form. In doing this, it is necessary to introduce canonical conjugate momentum densities 
which are defined by [9,10,14] 

^' = = ^7° = ^+ 
IT— — 9C — 

^ = IF = ~ dtA± ' (52) 
n = (jgfek = 

where the subscripts R and L mark the right and left-derivatives with respect to the real time respectively. With the 
above momentum densities, the Lagrangian in Eq. (51) can be represented as [9,10,14] 

c = n^v + wd t A^ + m t c + d t cu - n (53) 

where 

H = <M 7^ • v + m)v + i(n M ) 2 - ^„v% - nn + CV 2 C + iefa^Ap (54) 

is the Hamiltonian density. This Hamiltonian density is now written in the Euclidean metric for later convenience. 
The matrix 7^ in this metric is defined by 74 = [3 and 7 = — i(3a [14]. It should be noted that the conjugate quantities 
n and n for the ghost field are respectively defined by the right-derivative and the left one as shown in Eq. (52) 
because only in this way one can get correct results. This unusual definition originates from the peculiar property of 
the ghost fields which iare scalar fields, but subject to the commutation rule of Grassmann algebra. 

In order to derive the coherent-state representation of the action of thermal QED, one should employ the Fourier 
transformations for the canonical variables of QED which are listed below. For the fermion field [5,6, 10], 



^ r) = J [u s (p)b s (p, r)e^ + V s {p)d* s [p, rje"^ 2 ] , 

4(x , t)=J [u s (p)b* (p, T )e-* s + v s (p)d s (p, r)e^ 



where u s (p) and v s (p) are the spinor wave functions satisfying the normalization conditions u s+ (p)u s (p) = 
v s+ (p)v s (p) = 1, b s (p, t) and b*(p, r) are the eigenvalues of the fermion annihilation and creation operators b s (p, r) and 
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bf(p, t) which arc defined in the Heisenberg picture, d s (p, r) and d*(p, r) are the corresponding ones for antifermions. 
For the photon field [5,6,10], 

MZt) = / (-2% 7=1 A {k)[ax{k,ry % - s + a\{%,r)e- i% -% 

U,(x , r)=ij (^W^ ^(*)M*> r)e^ - a* x (k, r)^] 

where s^(le) is the polarization vector. The expression of n M (x, r) follows from the definition in Eq. (52) and is 
consistent with the Fourier representation of the free field. For the ghost fields which are always free for QED, we 
have 

r) = J [c{q, r)e^ S + c* (q, T )e-**\ , 

C(x, r) = J (2 $ /2 [c(q, T )e«* + c* (q, r)e~^] , ^ 



H(f , r)=ij ^^/f[c(J T)e** - c*(q, r)e-*% 



The expressions of 11(2?, r) and II(x, r) also follow from the definitions in Eq. (52). 

For simplifying the expressions of the Hamiltonian and action of the thermal QED, it is convenient to use abbrevi- 
ation notations. Define 

"'Mr) = J1:£:!' '/A = H (58) 



^H< 2 :»:^s-y/= + 7 <»> 



and furthermore, set a = (p,s,8) and 



then , Eq. (55) may be represented as 



d* s (p,r), if e 

)" 3 /V(p), if 
(2n)-3/ 2 v s (p), if 6 : 

£ = £/^ 



^(x,T)=EW a b* a (r)e- ie P- s . 



(60) 



(61) 



Similarly, when we define 



ja x (k,T), if 6 = +, 
a X e{k,T) = \ r \ (62) 



a* x (k,T), if 



n^(fc) = ^(2 7 r)- 3 / 2 [o;( g -)/2] 1 /2 e A ( fc ) 



and furthermore, set a = (k, A, 6*) and 



(63) 



(64) 



Eq. (56) can be written as 

A„(*,T) = £^aa(r)e* E -*, 
n„(^r) = En«a a (T)e < «K 



(65) 
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For the ghost field, if we define 



and furthermore set a = (q,a,9) and 



Eq. (57) will be expressed as 



9,~ \ J c a(<7», if 9 = +,\ 



E = E / d ^ ( 6 8) 



-iOq-x 



n»(f,T) = En a c C( (T)e i 
rf(f,r) = En a c;(r)e- 



(69) 



Now we are ready to derive the expression of the action given by the Lagrangian (53) in the coherent-state repre- 
sentation. For this purpose, it is convenient to use the following boundary conditions of the fields [5,6]: 



ip(x,0) = ip(x), V>(x,0) = ip(x), 
i>(x,0) = -ip{x), i>(x,/3) = -tp(x), 



(70) 



A„(x,0) =A ll (x,/3) =A„(f), , 
n„(j?,0) = !!„(:?,/?) =n„(f) {n > 



and 



C_(x, 0) = C(x, (3) = C(x), C(x, 0) = C(x, (3) = C(x), 
IL(x, 0) = U(x, (3) = U(x), U(x, 0) = IL(x, (3) = IL(x). 

Based on these boundary conditions, by partial integrations, the action may be written in the form 

+| [H„ [x, t)A^ (x, t) - n„ {x, t)A^ (x, t)] 
+f[II(x, r)(7(f, r) - n(f, t)C(x, t) 

+C{x, t) II (x, t)- C 7 (x, t)TT(x, t)] - H{x, t)} 



(72) 



(73) 



where the symbol " • " in ip(x, r), A^(x, r) • --denotes the derivative of the fields with respect to the imaginary time r. 
Upon substituting Eqs. (61), (65) and (69) into Eq. (73), it is not difficult to get 

S = - Jo dT U d 3 H±[b*(k, T%(k, t) - b* s (k, T)b s 0, t)} + ±[d*(k, T)d s (k, t) 

-d* s (k,T)d s (k,T)] + \ [a\{k,T)a x (k,T) - a* x (k,T)a x (k,T)} + ±[c*(fc,r) c (fc,r) ^ 

- c* (k, r)c{k, t) - c* (k, r)c(fc, t) + c* (k, r)c(k, r)]} + H(t)} 

= —Se 

where Se is the action defined in the Euclidean metric and 

H(t) = J d 3 xH(x) = H (t) + i?j(r) (75) 

in which Hq(t) and Hi{t) are respectively the free and interaction Hamiltonians given in the coherent state repre- 
sentation. Their expressions are 
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H (r)=j:O a e a b*jT)b a (T) 
+lEw«<WMT) + EwaC*(r)c a (T) (76) 

a a. 

and 

H i( T ) = E /("^^aW^W^W (77) 

a/37 

where 

/(a07) = ie(2ir) 3 6 3 (6 a p a - 6pp - O^W^ {p a )l^ (p^A 6 ^ (fc 7 ). (78) 

In the above, 8 a = 6, e a = (p 2 + M 2 ) 1 / 2 is the fermion energy, to a = k is the energy for a photon or a ghost particle. 

It is emphasized that the expressions in Eqs. (75)-(78) are just the Hamiltonian of QED appearing in the path-integral 
shown in Eq. (42) where all the creation and annihilation operators in the Hamiltonian which are written in a normal 
product are replaced by their eigenvalues. It is noted that if one considers a grand canonical ensemble of QED, the 
Hamiltonian in Eq. (74) should be replaced by K(t) defined in Eq. (2). Employing the abbreviation notation as 
denoted in Eqs. (58), (62) and (66) and letting q a stand for (a a ,b a , c a ), the action may compactly be represented as 

S E = f dT{J2l[Q* a (r)oq a (r)-q* a oq a (T)}+H(r)} (79) 
Jo a * 

where we have defined 

q* a °q a = a a - a a + + b* a b a + 9 a c* a c a (80) 

and the Hamiltonian was represented in Eqs. (75)-(78). It is emphasized that the a = ± is now contained in the 
subscript a. Therefore, each a may takes a + and/or a" as the first term in Eq. (80) does. 
Now, let us turn to the <p 4 theory whose Lagrangian at zero-temperature is 

C = \d»<pd& - imV - (81) 

With the canonical conjugate momentum density defined by 

dC 

H = _ = I* = vp t (82) 

the Lagrangian in Eq. (81) at finite temperature may be written as 

C =nd t ip - H (83) 

where 

H = \^ + ^(V,) 2 + ImV + JV- (84) 

To derive the coherent-state representation of the action given by the Lagrangian in Eq. (83), we need the following 
Fourier representation of the canonical variables: 

<P(2, r) = J -7=- Hk, T)e&* + a* (k, r)e~^\ , 

n(x, r) = -if yj Sfl [a(k, r)e ik -* - a* (k, r)e~^\ 

where to(k) = (k 2 + m 2 ) 1 / 2 . When we define 

* ( *' T H«-(*V), = (86) 
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9e(k) 



5 e (fc) = H) (5^\/# 



(87) 



and furthermore set a = (k, 6) and 



Eq. (85) can be written as 



E=/ d3fc E> ( 88 ) 



With these expressions and considering the boundary conditions 



ir(x, (3) — ir(x, 0) = n(x), 



the action given by the Lagrangian in Eq. (83) can be found to be 

S = J Q dr J d 3 x{^[n(x, r)<p(x, t) — -k{x, r)ip(x, t)] — H(x, t)} 



where 



in which 



So dT Ul \ K (r)a a (r) - d* (r)a a (r)] + H(r)} (91) 

a 

= — Se 

H(t) = H q (t) + H t (t) (92) 

H (t) =E w « a --( T ) a «+( T ) ( 93 ) 



is the free Hamiltonian and 



H i( T ) = E (2n) 3 5 3 (6 a k a + Opkp + 9 y k 7 + 6sks)^g a gpg-yg8a a a/3a 7 as (94) 
is the interaction Hamiltonian. 

IV. GENERATING FUNCTIONAL OF GREEN'S FUNCTIONS 

With the action Se given in the preceding section, the quantization of the thermal fields in the coherent-state 
representation is easily implemented by writing out its generating functional of thermal Green's functions. According 
to the general formula shown in Eq. (42), the QED generating functional can be formulated as 

ZV]=fDtfq)e-*'-*fV(q*q)exp{h[f(0)-Q(0) m] 
- q *(0)- q (0)}-SE + Sodrr(T)- q (r)} [ > 



where we have defined 



and 



q* ■ q = ^a* a a a + 9 a b* a b a + c* a c a (96) 



3* ■ 1 = C a a a + e a ( V * a b a + b* aVa + Qc a + c* a ( a ) (97) 
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here £, a ,Va an d Ca are the sources for photon, fermion and ghost particle respectively and the repeated index implies 
summation. It is noted that the product q* ■ q defined above is different from the q* a o q a defined in Eq. (80) in the 
terms for fermion and ghost particle and the subscript a in Eqs. (96) and (97) is also defined by containing 9 a = ± 
. In what follows, we assign to represent the a with 6 a — ±. According to this notation, the sources in Eq. (97) 
are specifically defined as follows: 

£* 

V* a , (98) 

Sa 

where the subscript a on the right hand side of each equality no longer contains a and the photon term in Eq. (96) 
(l/2)a*a Q may be replaced by a a -a a +. The integration measures D{q*q) and D(q*q) arc defined as shown in Eqs. 
(24) and (25).' 

Now we are interested in describing the perturbation method of calculating the generating functional. Since the 
Hamiltonian can be split into two parts H (t) and H[(t) as shown in Eqs. (76) and (77), the generating functional 
in Eq. (95) may be perturbatively represented in the form 

Z\j] = cxp{- J P drff 7 (^y)}Z°[j] (99) 

where Z°[j] is the generating functional for the free system and the exponential functional may be expanded in a 
Taylor series. In the above, the commutativity between Hi and Z°[j] has been considered. Obviously, the Z°[j] can 
be written as 



Sq+ Sen Sc« 
Va+ = Va, Va- = 
Ca+ — Ca, Cct~ = 



Z°\j] = Z»[(]Z»[rj\Z%] (100) 

where Zp[£], Z®\rj\ and Z®[Q are the generating functional contributed from the free Hamiltonians of photons, 
fcrmions and ghost particles respectively They are separately and specifically described below. 

In view of the expressions in Eqs. (95), (74) and (76), the generating functional is of the form 



= I D ( a * a ) ex P{- / d 3 k[a* x (k)a x (k)} 
x jD(a*a) exp{I p (a* x ,ax-,tx,£\)} 



(101) 



where 



h to , a* ; & ) = / d 3 k\ [a* x (k, P)ax (£,/?) + a* x (k, 0)a x (k, 0)] 
dr J d 3 fc{i[a A (fc, r)a x (k, t) - d* x (k, r)a x (k, r)] + w{k)a\(k, r)a x (k, r) (102) 
-Cx(k,r)ax(k,r)-a* x (k,r^x(k,r)} 



and 



D(a*a) = l\±da* x (k)dax(k), 



kX 



®(«*«) - II ±da* x (k,T)dax(k,T). 



(103) 



fcAT 



The subscript A in the above denotes the polarization. When we perform a partial integration, Eq. (102) becomes 

Ip{a* x ,a\;Z\,Z\) = J d 3 ka* x {k,f3)a x {k,f3) - j£ dr J d 3 k{a* x {k,T)d x {k,T) 
+ W (fc)a A (fc,r)a A (fc,T)-^(fc,r)a A (fc,r)-a A (fc,T)6(^r)}. 

For the generating functional Z° [rj\ , we can write 

z° f [ V ] = jD(b*bd*d)cM-fd 3 k[b:(k)b s (k)+d*(k)d s m (105) 

x / D (b* bd* d) exp{// (6* , b s , d* s , d s ; rf s , iq s , rf s ,rj„ ) } 

where 
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I f (b*„ b s , d* s , d s ■ V * s , Vs , rj* s ,Vs ) = S d 3 k \ IK (k, 0) b s (k,0) + d* s (k,0)d s (k, 0) 
+K(k,0)b s (k,0) + d*(k,0)d s (k,0)}- dr J d 3 k{±[b* s (k,T)b s (k,r) 

-b* s (k, r)b s (k, r)] + \ [d* s (k, r)d s (k, t) - d* s (k, r)d s (k, r)] (106) 
+e(k) [b* s (k, r)b s (k, t) + d* s (k, r)d s (k, r)] - [ V * s (k, r)b s (k, r) 
+b* s (k, t) Vs (k, t) + rj* s (k, r)d s (k, t) + d* s (k, T )rj s k, r)]} 



and 



D(b*bd*d) = l\db* s (k)db s (k)dd* s {k)dd s (k), 

£s -. -. - - (107) 

D(b*bd*d) = Hdb* s (k,T)db s (k,T)dd* s (k,T)dds(k,T) v ' 

ksr 

in which the subscript s stands for the fermion spin. By a partial integration over r, Eq. (106) may be given a simpler 
expression 

If (b* s ,b s ,d* s ,d s ;r,* s ,Vs,K,rj s ) = J d 3 k[b* s (k, 0)b s (k,0) + d* s (k, 0)d s (k, 0)} 
- jj dr J d 3 k{b* (k, t% (k, t) + d* s (k, r)d s (k, t) 
+e(k) [b* s (k, r)b s (k, t) + d* (k, r)d s (k, r)} - [77* (k, r)b s (k, r) 
+b* s (k, t) Vs (k, t) + % (k, r)d s (k, t) + d* s (k, t)tJ s k, t)}}. 

As for the generating functional Zg [£] , we have 

Z° g [(} = J D(c*ccc*)cxp{- J d 3 k[c*(k)c(k) - c*(k)c(k)]} 
x J* J) {c*ccc* ) cxp{/ 9 (c* , c, c* , c; C* , C, C* , 0} 



(109) 



where 



I g (c* , c, c* , c; C* , C, C* , = / <^5 [c* (fc, P)c(k, 0) - c* (k, 0)c(k, 0) 
+c* (jfe, 0)c(k, 0) - c* (jfe, 0)c(fc, 0)] - /„" dr / d 3 fc{ ± [c* (jfe, r) c (jfe, r) 

-c* (fc,r)e(fc,T)]-|[c*(jfe,T)c(jfe,r)-c*(jfe,T)c(ife,r)] ( 110 ) 
+w(fc)[c*(fc,r)c(fc, r)_-c*(A,r)c(A, r)] - [C*(ft,r)c(ft, r) 
+c*(fc, r)C(fc, r) + (*(k, r)c(k, t) + c* (k, r)C(A, r)]} 



and 



D(c*ccc*) = n^*(fc)<^(^)rfc(fc)dc*(fc), 
S(c*ccc*) = [J dc* {k, T)dc(k, r)dc{k, T)dc* (jfe, t). 

feT 



(111) 



After a partial integration, Eq. (110) is reduced to 



I g (c* , c, c* , c; C* , C, C*, = / d 3 k[c* (k, 0)c(k, 0) - c* (k, 0)c(k, 0) 
— J J d 3 k{c*(k, r) c (fc, r) — c*(k, r)c(fc, t) + w(fc)[c*(fc, r)c(fc, r) 
-c* (fc, r)c(fc, r)] - [C* (fc, r)c(fc, r)_+ c* (fc, r)C(fc, r) + f (*, r)c(fe, r) 
+c*(fc,r)C(fc,r)]}. 



(112) 



Here it is noted that all the terms related to the quantities c* and c are opposite in sign to the terms related to the 
c* and c and, correspondingly, the definitions of the integration measures for these quantities, as shown in Eq. (Ill), 
are different from each other in the order of the differentials. 

The generating functionals in Eqs. (101), (105) and (109) are of Gaussian-type, therefore, they can exactly be 
calculated by the stationary-phase method. First, we calculate the functional integral Zp[£\. According to the 
stationary-phase method, the functional Zp[£] can be represented in the form 

Z° p m - J D(a*a)exp{-J d 3 k[a* x (k)a x (k) + I°(a* x , a x - ft, 6)]} (113) 



13 



where 7° (a x , a x ; £ x , £ A ) is given by the stationary condition <57 p (a A , a A ; £\) = 0. By this condition and considering 
the boundary conditions [4-6]: 

a* x (k,(3) = a* x (k),a x (k,0)=a x (k), (114) 
one may derive from Eq. (102) or (104) the following inhomogeneous equations of motion [7,9]: 

a x (k, t) + u(k)a x (k, t) = r), . . 

d* x (k,r) - u;(k)al(k,T) = -C x (Kt). 1 ' 

In accordance with the general method of solving such a kind of equations, one may first solve the homogeneous linear 
equations as written in Eq. (33). Based on the solutions shown in Eq. (34) and the boundary condition denoted in 
Eq. (114), one may assume [7,9] 

a x (k, t) = [a x (k) + u x {k, r)]e~"( fc >, 
al(k,T) = [al(k)+u* x (k,r)}e^Kr-P) 

where the unknown functions u x (k,r) and u* x (k,r) are required to satisfy the boundary conditions [7,9]: 

u x (k, 0) = u x (k, /3) = u* x (k, 0) = u* x (k, (3) = 0. (117) 
Inserting Eq. (116) into Eq. (115), we find 



«a(t) =&(fc,T)e<"(*> T ) 
u* x (t) = -C x (k,T)e"WW- T \ 

Integrating these two equations and applying the boundary conditions in Eq. (117), one can get 

u x (k,r)= f dT'e^y^{ky), 

<(k,r) = -j;dr'e^-^C X (k,r'). 
Substitution of these solutions into Eq. (116) yields [7,9] 

a x (k,r) = a x (k)e-^ + £ dT'e-^^-^Uk,r'), 
a* x (k,r) = al(k)e^ T -V + J? dT'e^ T - T '^l(k,T'). 



(118) 



(119) 



(120) 



When Eq. (120) is inserted into Eq. (102) or Eq. (104), one may obtain the 7° (a* x , a x : £ x , £\) which leads to an 
expression of Eq. (113) such that 

z p[i\ = J D(a* a) cxp{- J d 3 k[a* x (k)a x (k)(l - e~^) - a* x {k)e~^ 

x/ ^re-( fe >a(fc,r)-/ /3 dr e --( fe >a(fc,T) aA (fc)] (121) 
+ So dri So ^ / d 3 k? x (k, n)fl(ri - T2 )e-^)(rr-r T2 ) } . 

The above integral over a* x (k) and a x (k) can easily be calculated by applying the integration formula denoted in Eq. 
(37) when we set 

A = 1 _ e -/Mfe) 5 

b = e-P"® J o dre^ T ^ x {k, t), (122) 
f(a)=J Q P dre-^)r Cx{ k jT)ax{ k). 



The result is 



where 



Zl\i] = Z° p cxp{- f dn f dr 2 [ d 3 k^ x * (k, n)A AV (k, n - r 2 )£ A ' (k, r 2 )} (123) 
Jo Jo J 
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Z p = II t 1 - e-^}- 1 = - e -^}- 4 (124) 

which is precisely the partition function contributed from the free photons [5,6] and 

A AV (fc,Ti - r 2 ) = AV A 6 (fc,ri - r 2 ) (125) 

with 

A 6 (jfe,ri - r 2 ) = [fl(ri - r 2 ) - (1 - ^)-i] e -»(i)(^) (126) 

is the free photon propagator given in the Feynman gauge and in the Minkowski metric {Note: in Euclidean metric, 
g\\t — > — <Jaa')- When we interchange the integration variables n and t 2 and make a transformation k — > — fc in Eq. 
(123), by considering the relation 

r) (127) 

which will be interpreted in the next section, one may find that the propagator in Eq. (126) can be represented in 
the form 

A b (k, Tl -T 2 ) = i[n fc (fc)e^ (£)|Tl - T21 -n 6 (fc)e" (£)|ri - T21 ] (128) 

where 

n b (k) = (1 - e-^)-\ n b {k) = (1 - e^)" 1 (129) 

which are just the boson distribution functions [3-6,13]. 

Let us turn to calculation of the functional integral in Eq. (105). On the basis of stationary-phase method, we can 
write 

z° f [ v ] = J^(6*6d*d)cxp{-/d 3 fcK(fc)Mfc)_+d:(fc)4(fc)]} (130) 

x exp{I^(b*,b s ,d* s ,d s ;r]* s ,r] s ,r]* s ,r] s )} 

where If(b* s , b s , d* s , d s ; rf s , r} s ,rf* a , rj s ) will be obtained from Eq. (106) or Eq. (108) by the stationary condition 
SIf(b*,b a ,d*,d 8 ;T)*,rf a ,Tj*,Tj s ) = 0. From this condition and the boundary conditions [4-6] 

b* s (k,(3) = ~b* s (k), b.(k,0) = b s (k), 

d* s (k,p) = -d* s (k), d s (k,0) = d s (k), [ ' 



one may deduce from Eq. (106) or Eq. (108) the following equations [7,9] 

(132) 



b s (k,r) +e(k)b s (k 1 T) = r) s (k,j), 
b\{k,j) - e(k)_b*(k,j) = -rj*(k,T), 

d s (k,r) + e(k)d s (k,r) = rj 3 (k,T), 
dt(k,r)-6(k)d* s (k,T) = -r}* s (k,r). 



Following the procedure described in Eqs. (115)-(120), the solutions to the above equations, which satisfies the 
boundary conditions in Eq. (131) and the conditions like those in Eq. (117), can be found to be [7,9] 

b s (k,r) = b s {k)e- £ ^ T + J* r dT'e-^^-^rjsik,^), 
b* s (k,r) = -b* s (k)e^-V +^dr'e^- '^(fc, r'), , , 

d s (k, t) = d s {k)e- e ^ T + J T dr'e- £ ^ r - T '\{k, r'), 
d*(k,r) = -d* s {k)e e ^ iT -^ + J? dT'e^ T - T ")r]* s (k,T'). 

Substituting the above solutions into Eq. (106) or (108), we find 

J° (b* s ,b s ,d* s ,d s ;r]* s , Vs ,rl* s ,r] s ) = Jd 3 k{-e~^) [b* (k)b s (k) + d* (k)d s (k)} 

+ / /3 dre-^[^(fc,T)6 s (fc) +rj*s(k,r)d s {k)} - J Q dre^ T (134) 

x [b* (k) Vs (k,r)+ d* (k)rj. (k,r)]} + B [rf t , Vs , r£ , rj s ] 
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where 



B[ V *, Vs ,r s ,r} s ] = \ If dn dr 2 J d 3 k{9( Tl T2 ) e -e(k)(r,-r 2 ) 



x[r)* s (k,T 1 )r ls (k,T2) +r s (k J T 1 )rj s (k,r 2 )} + 0(r 2 - Tl ) e -c(*0(— (135) 
x [77* (fc, T 2 )Vs(k,n) + rj* s (fc, t 2 )rj s {k, n )] }. 

On inserting Eq. (134) into Eq. (130), we have 

Z° f [ V ] = A^^sX.Vs]^^ 111 ^ (136) 

where 

A[ V * s ,VsJfs,ris} = I D(b*b) cxp{- / d 3 k[b* s (k)b s (k)(l + e-^(£)) 

x / £>(<i*d) exp{- / d 3 fc[^(fc)d s (fc)(l + e-^( fe )) + e -^«d*(fc) 
x dre<^rj s (k, r) dre"'^ (fc, r)d s (k)}} 

here the fact that the two integrals over {6*, b} and {d* , d} commute with each other has been considered. Obviously, 
each of the above integrals can easily be calculated by applying the integration formulas shown in Eq. (39). The 
result is 

A[ V t, Vs X s ,Vs] = ZjeM-Uodnjy^Jd^l + e^)^ 

x {e~<^-^ [ v *(k, Tl ) Va (k, t 2 ) + t s (k, n)Tl s {k, r 2 )\ (138) 

+e -e(k)(r 2 -r l)[vt ( £ ^fafc n ) + r 2 )r) s {k, Tl )}}} 

where 

Z° f = J][l + e-^)] 2 (139) 

which just is the partition function contributed from free fermions and antifermions [4-6,13]. It is noted that the two 
terms in the exponent are equal to one another as seen from the interchange of the integration variables t\ and t 2 . 
After Eqs. (135) and (138) are substituted in Eq. (136), we get 

Zj[rj\ = Z° f exp{jj dn dr 2 J d 3 k{[9( Tl - r 2 ) - (1 + e ^W)-i] e -e(fe)(r 1 -r 2 ) 

xlnKk^Mk,^) + rf s {k,T 1 )l) s {k,T 2 )] + [6{t 2 - n) - (1 + e^)- 1 } ( 140 ) 

xe -e(fc)(r 2 -r l)[r? * ( fc ; ^(^ n) + ^ ^(fc, n) ]} 

When we interchange the variables n and r 2 , set fc — > — fcin the second term of the above integrals and note the 
relation 

V*s(k,T2)Vs(k,Ti) = r?*(-fc,ri)r? 5 (-fc,T 2 ) (141) 
which will be proved in the next section, the functional Z®[rj\ will eventually be represented as 

Z° f [ry] = Z° f eMIo dT i Jo d ^ I d3k ^t (k, n) AJ (fc, n - t 2 )v (fc, r 2 ) 



where 



ith 



Af (fc, ri - r 2 ) = 5 SS A/(fc, ri - r 2 ) (143) 



A/^n - t 2 ) = i[n/(fc)e- £ ( £ )l Tl - T2 l - ri / (fcV (£)|Tl ~ r21 ] (!44) 



being the free fcrmion propagator in which 
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n f (k) = (1 + e"^)- 1 , n f (k) = (1 + e^)" 1 (145) 

are the familiar fermion distribution functions [4-6, 13]. 

Finally, let us calculate the generating functional Z g [£] . From the stationary condition 51 g (c* ,c,c* , c; £* , C> C > C) = 
and the boundary conditions: 

c*(k,f3) = c*(k), c*(k,p) = c*(k),c(k,0) = c(k),c(k,0) = c(k) (146) 

which are the same as those for scalar fields other than for fermion fields [5], it is easy to derive from Eq. (110) or 
(112) the following equations of motion 

c(k,r)+w(k)c{k,r) = -C(£,t), 
c*(fc,r) - uj(k)c* a (k,T) = Q(k,r), 

c{k,r)+u{k)c a {k,r) =C„(fc,r), { ^' } 
c* (k,r)-oj(k)cl(k,T) = -(* a (k,T). 

By the same procedure as stated in Eqs. (115)-(120), the solutions to the above equations can be found to be 



(148) 



c(k, t) = c(k)e-^ T - J T dT'e-^)( T - T \(fc, r'), 
c*{k,r) = c*(k)e^ T -^ - J? dT'e^ T - T '\*{k,T'), 

c(k, t) = c{k)e-^ T + J r dT'e-^^'^ik, r'), 
c*(k,r) = c*(k)e^ T -^ + J? dT'e^ T - T '^*{k,T'). 

Upon substituting the above solutions into Eq. (110) or Eq. (112), we find 

I° g (c* , c, c* , c; C, C, C* , C) = / d 3 fc{e-^( fe ~) [c* (k)c(k) - c* (fc)c(fc)] 

+e -^a»(S) dTe "(fe)r [ c * (fc)C(fc, r) + C* (fc)C(fc, T)] 

- So dr i So dT ? S d3k0 (ji - t^c^-t.) [c Ti)c( ^ T2) 

-C*(fc,Tl)C(£,T 2 )]. 

On inserting the above expression into the following integral given by the stationary-phase method: 

Z° g [(} = J D(c*ccc*)eM- J d 3 k[r(knk)-c*(k)c(k)} + 7°(cV, c*, c; C*, C,f , 0} (150) 
and applying the integration formulas in Eq. (39), one can get 

Z° g [Q = ZVeMSodnSo^Sd'mn r 2 ) - (1 - e /M£))-i] e -«(£)(n-T*) 

x[c*(fc,r 1 )C(fc,r 2 )-C(fc,n)C(fc,T 2 )]} 

where 

Z°==n[l-e^ (£) ] 2 (152) 
t 

is just the partition function arising from the free ghost particles which plays the role of cancelling out the unphysical 
contribution contained in Eq. (124). If we change the integration variables in Eq. (151) and considering the relations 

C(k,r) = -C(-k,r), ((k,r) = -f(-fc,r) (153) 
which will be interpreted in the next section, Eq. (151) may be recast in the form 

Z° g [(] = Z° g eMSo dT ^ So dr 2 / d 3 k[(*(k, n)A b (fc, n - r 2 )C(k, r 2 ) 
-Cnfc,n)A 6 (fc,Ti-T 2 )C(fc,T 2 )]} 

where Ab(k,n — r 2 ) was written in Eq. (128). 
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Up to the present, the perturbative expansion of the thermal QED generating functional in the coherent-state 
representation has exactly been obtained by the combination of Eqs. (99), (100), (123), (142) and (154). Especially, 
the partition function for the free system has been given by the product of Eqs. (124), (139) and (152). The partition 
function for the interacting system can be calculated in the way as shown in Eq. (45). Here it should be noted that 
the differential 6/8 j(r) in Eq. (99) represents the collection of the differentials 8/8^(k,r), 8/S^\(k,r), —S/8r] s (k,T), 
5/Sr,t(k,T), -8/8r) s {k,r), 8/8rj*{k,r), -8/8((k,T), 8/SC(k,r), -8/8((k,T) and 5/5?(k,T). 

In the last part of this section, we briefly describe the generating functional for the thermal if 4 theory. The generating 
functional may be written out from Eq. (95) by letting q = a. Its perturbative expansion is still represented by Eq. 

I 5 > 
v<y(r)> 



(99) with noting that the H I (j^- 1 ) is given by Eq. (94) with replacing a a in Eq. (94) by S/Sj* and 



Z [j] = J D(a*a) cxp{- / d 3 ka* (k)a(k)} J®(a*a) exp{/(a*, a;j*,j)]} 
= J D(a*a) exp{— J d rka* (k)a(k) + Io(a* , a; j* , j)} 



where 

r/3 



I {a* , a; j* ,j) = J d 3 ka* (k, 0)a(k, [3) - ^ dr J d 3 k{a* (fc, r)d(fc, r) 
+cu(k)a*(k, r)a(fc,r) -j*(k,r)a(k, r) - a*(k, r)j(k, r)} 



(155) 



(156) 



and Io(a* , a; j* , j) is determined by the stationary condition 8I(a*,a;j*,j) = 0. Completely following the procedure 
described in Eqs. (113)-(129) with noting the boundary condition 

a* (k, (3) = a (k), a(k, 0) = a(k) (157) 

and the relation 

j*(k,r)=j(-k,r), (158) 
which will be proved in the next section, it is easy to obtain 

ZoK]=^°exp{/ dn dr 2 d 3 kf(k,T 1 )A b (k,T 1 -T 2 )j(k,T 2 )} (159) 
Jo Jo J 



where 



7 o 



= TT (i6o) 

k 

and Ab(k, t\ — t 2 ) was represented in Eq. (128). 



V. DERIVATION OF THE GENERATING FUNCTIONAL REPRESENTED IN THE POSITION SPACE 

In this section, we plan to show how the perturbative expansions of the generating functionals represented in the 
position space are derived from the corresponding ones given in the preceding section. For this purpose, we need to 
derive the generating functional represented in position space for free systems. For thermal QED, it can be written 
as 

Z°[J] - Z° f [I,I]Z^]Z° g [K,K] (161) 

where Z® [I, I] , Z° [ J M ] and Z° [K, K] are the position space generating functionals arising respectively from the free 
fermions, photons and ghost particles and /, /, J M , K and K are the sources coupled to the fermion, photon and ghost 
particle fields, respectively. In order to write out the Z*j [1,1], Zp[J/j] and Zg[K,K] from the generating functionals 
given in Eqs. (142), (123) and (154), it is necessary to establish relations between the sources introduced in the position 
space and the ones in the coherent-state representation. Let us separately discuss the functionals Z® [1,1], Z^J^] and 
Zg[K,K]. First we focus our attention on the functioinal Z®[I,I]. Usually, the external source terms of fermions in 

the generating functional given in the position space are of the form J Q dr J d 3 x[I(x, r)ij)(x, r) + ip(x, t)I(x, t)] [5,6]. 
Substituting in this form the Fourier expansions in Eq. (55) for the fermion field and in the following for the sources 
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H^r) = J T A^I(k,r)e ik -* 
7(Z,T)=fj£^I(k,T)e-^ 



(162) 



we have 



where 



dr J d 3 x[I(x, t)1>{2, t) + <P(x , t)I(x, t)} = f£ dr J d 3 k[ V * s (k, r)b s (k, r) 
+b* s (k, T )vs (k,r)+ fj* s (k, r)d s (k,r) + d* s (k, r)rj s (k,r)} 



r) s (,r) = u s {k)I{k,T), rf s {k,r) = I{k,r)u s {k), . . 

r] s (k,r) = -I(-k,r)v s (k),r s (k,r) = -v s (k)I(-k,T). [ ' 

From these relations and the property of Dirac spinors, the relation in Eq. (141) is easily obtained. On substituting 
Eq. (164) into Eq. (142), one can get 

Z° f [I, I] = Z) cxp{ / dn / dr 2 / d 3 kl{k, n)S F (k, n - r 2 )I{k, t 2 )} (165) 
Jo Jo J 

where 

S F {k,n-T2) = [Qt + M)/s(k)}A f (k,T 1 -T 2 ) (166) 

with k =7 jU fc'' and fc M = (k, ie n ). By making use of the inverse transformation of Eq. (162), the generating functional 
in Eq. (165) is finally represented as [5,6] 

r0 rP 

Z° f [I,I] = Z° f exp{ d 4 Xl d i x 2 I(x 1 )S F (x 1 - x 2 )I(x 2 )} (167) 
Jo Jo 

where x = (x, r), d 4 x = drd 3 x and 

C d 3 k 

S F ( Xl -x 2 )= / -^S F (k,n - r 2 )e lk - 3 (168) 



, (2tt) 3 

It is well-known that the propagator Af(k,T\ — t 2 ) is antiperiodic, 

A f (k, n - t 2 ) = -A f (k, n-r 2 - 0), if n y t 2 , lgg , 
A f (k,T 1 -T 2 ) = -A f (k,T 1 -T 2 + p), ifn <t 2 

This can easily be proved from its representation in the operator formalism shown in Eq. (46) with the help of 

the translation transformation b s (r) = e l3K b s e~ l3K and the trace property Tr(AB) = Tr(BA). According to the 
antiperiodic property of the propagator, we have the following expansion 

^f{k,r) = W^f{le n )e~^ (170) 

n 

where r — ri — T2, e n — j§(2n + 1) with n being the integer and 

A f (k, e n ) = f dre^A f (k, r) = ^ (171) 
Jo ei + eik) 2 

where the expression in Eq. (144) has been used. Substituting Eqs. (170) into Eq. (166) and noticing the above 
expression, the propagator in Eq. (168) can be written as [5,6] 

I r d 3 k £ i k ■(x 1 -X2)-ie n (r 1 - n2 ) 

S F {x! - x 2 = - > / —- = 172 

n J 2?r 7-k + M-ie n7 <> 
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Next, we discuss the generating functional Z^[J^]. The source term of photons in the generating functional given 

in the position space is commonly taken as dr J d 3 xJ^(x, t)A^(x, r) [5,6]. Employing the expansions in Eq. (56) 
and in the following 

/d 3 k 
^-^J^k,r)e ik - S , (173) 



we can write 



J*dr J d 3 xJ^x,r)A^x,r) = J*dT J d 3 k[C x (k,T)a x (k,T) + a* x (k,T)Uk,r)} (174) 

where 

6(fc,r) = (2u;(k))- 1 / 2 el(k)J,(k,T) - 8(-fc,r) (175) 

in which the last equality follows from that the J M (x, t) is a real function. Inserting the relation in Eq. (175) and 
then the inverse transformation of Eq. (173) into Eq. (123) and considering completeness of the polarization vectors, 
one may find the generating functional Z^J^] such that [5,6] 



1 f 13 t 

Z°[J^} = Z p °cxp{-- jf d A Xl J d 4 x 2 J^ Xl )D^{ Xl - x 2 )r(x 2 )} 



(176) 



where 

/d 3 k 1 
A 6 (fc,n -r 2 )e*^-^) (177) 
oj{k) 

By the same argument as mentioned for Eq. (169), it can be proved that the photon propagator A(,(fc, n — t 2 ) is a 
periodic function 

A b (fc,Ti - r 2 ) = A(,(fc,n - r 2 - /?), i/ ri >- r 2 ; ^ 17g ^ 
A b (fc, ti - r 2 ) = A 6 (fe, n - t 2 + /?), i/ n -< t 2 . 

Therefore, we have the expansion 

= 4 E A ^ w «) e ^" T ( l79 ) 

' n 

where co n = and 

A h (k, u n ) = f dre^A b (k, r) = - ^ (180) 

where the expression in Eq. (128) has been employed. Upon substituting Eqs. (179) and (180) in Eq. (177), we 
arrive at [5,6] 

D fn To) - 1 V f rf3fc ^ ifc-(xi-X2)-ia, n (Ti-T 2 ) QgjN 

M " ( 1 ''P^J {^) 3 el+.{kY e (181) 

which just is the photon propagator given in the position space and in the Feynman gauge. 

Finally, we turn to the generating functional Z®[K, K]. In accordance with the expansions in Eq. (57) for the ghost 
particle fields and those for the external sources: 



K(x,r) = J ^K(k,r)e ik - S 



tk-x 



(182) 



the relation between the sources in the position space and the ones in the coherent-state representation can be found 
to be 
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dr f <Px\K(2,t)C(2,t) + C(x,t)K(x,t)} = dr J_d 3 k[C(k,T)c(k,r) , . 

+ C*(k,T)C(k,T)+C*(k,T)c(k,T)+C*(k,T)C(k,T)} 

where 

C(fc,r) = (2u,(k))-i/*K(k,r), C(k,r) = (2u>(k))-^K(k,r), 
C(k,T) = -(2cj(kj)-^K(-k,r) 7 C(k,r)=-{Mk))- 1/2 K(-k,T). K ] 

from which the relations denoted in Eq. (153) can directly be deduced. When the above relations and the inverse 
transformations of Eq. (182) are inserted into Eq. (154), one can get 

Z° g [K,K] = Z° g cxp{- [ d A Xl [ d i x 2 K(x 1 )A g {x 1 - x 2 )K{x 2 )} (185) 
Jo Jo 



where 

AJxi -xo) = fA- , . 

which just is the free ghost particle propagator given in the position space [5,6]. In the derivation of the last equality 
of Eq. (186), the expansion given in Eqs. (179) and (180) have been used. 

With the generating functionals given in Eqs. (167), (176) and (185), the zeroth-order generating functional in 
Eq. (161) is explicitly represented in terms of the propagators and external sources. Clearly, the exact generating 
functional can immediately be written out from Eq. (99) as shown in the following 

Z[J] = exp{- jf ^ W/( _^_ )}z0[J] (lg7) 

where J stands for I, 7, J u K and K, represents the differentials - ffe, jj^-y and —gjfe and 

T~ti( s /( x ) ) can be written out from the last term in Eq. (54) when the field functions in the term are replaced by the 
differentials with respect to the corresponding sources. 

At last, we sketch the derivation of the generating functional given in the position space for the thermal ip theory. 
By using the expansions in Eq. (85) for the scalar field and in the following for the source, 

/ii 3 k 
7^)372 J(k,r)e ik -*, (188) 

it is easy to find the relation between the both sources given in the position space and the coherent-state representation 

r t r - 

dr d 3 xJ(x,T)ip(x,T) = / dr d 3 k[j*(k,T)a(k,T) +a*(k,T)j(k,T)} (189) 



where 

j(k, t) = {2uj(k))- 1 ' 2 J{k, t) = (2u(k)r 1 / 2 J* (-k, r) = f (-k, t) (190) 

The second equality in the above follows from the real character of the source J(x,t). When the above relation and 
the inverse transformation of Eq. (188) are inserted into Eq. (159), the generating functional Z g [J] is found to be 
[5,6] 

i r 13 r 13 

Z° g [J] =Z°exp{- J d 4 Xl J d i x 2 J(x 1 )A(x 1 -x 2 )J(x 2 )} (191) 



in which 

Mxi -x 2 ) = f ,fA ,V A b (fc, n - T2 ) e 'fc-(^-^) 

(192) 



A( Xl -x 2 ) = J ^-^A 6 (fc,n - r 2 )e ik <^-^ 

- I V f d3,C .. 1 _ r ik-(xi—X2) — iu) n (Tl— T?) 

(3 2— in J l2¥F e 2 +w( fe ) 2 t; 



where the expressions in Eqs. (179) and (180) have been considered. The exact generating functional is still represented 
by Eq. (187), but the Hi ( s f^ ) in Eq. (187) is now given by the last term in Eq. (84) with replacing ip(x) by jjj^- 
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VI. CONCLUDING REMARKS 

In this paper, the path-integral formalism of the thermal QED and (p 4 theory has been correctly established in 
the coherent-state representation. In contrast to the ordinary path-integral formalism set up in the position space, 
the expressions of the generating functionals presented in this paper not only gives an alternative quantization for 
these theories, but also provides an alternative and general method of calculating the partition functions, the thermal 
Green's functions and thereby other statistical quantities in the coherent-state representation. In particular, the 
generating functional enables us to carry out analytical calculations without concerning its original discretized form. 
The discretized form must be used for the generating functionals given in the previous literature because the previous 
generating functionals given in the coherent-state representation is incorrect and actually useless. As one has seen 
from Sect. 4, the analytical calculations of the zero-order generating functional is more simple and direct than the 
previous calculations performed in the discretized form either in the coherent-state representation or in the position 
space [4-6] . The coherent-state path- integral formalism established in the coherent-state representation corresponds to 
the operator formalism formulated in terms of creation and annihilation operators. In comparison with the operator 
formalism which was widely applied in the many-body theory [4,13], the coherent-state path- integral formalism has a 
prominent advantage that in the calculations within this formalism, use of the operator commutators and the Wick 
theorem is completely avoided. Therefore, it is more convenient for practical applications. It should be mentioned that 
although the QED generating functional is derived in the Feynman gauge, the result is exact. This is because QED 
is a gauge-independent theory. As shown in Sect. 4, in the partition function shown in Eq. (124) which is derived 
in the Feynman gauge, the contribution given by the unphysical degrees of freedom, i.e., the time and longitudinal 
polarizations of photons is completely cancelled out by the partition function shown in Eq. (152) which is contributed 
from the unphysical ghost particles. This indicates the necessity of retaining the ghost term in the effective Lagrangian. 
Certainly, the generating functional formulated in the coherent-state representation can be established in arbitrary 
gauges. But, in this case, the photon propagator would have a rather complicated form due to that the longitudinal 
part of the propagator will involve the photon polarization vector. Another point we would like to mention is that to 
formulate the quantization of the thermal QED and (p 4 theory in the coherent-state representation, we limit oursclf 
to work in the imaginary-time formalism. It is no doubt that the theory can equally be described in the real-time 
formalism. Needless to say, the coherent-state formalism described in this paper can readily be extended to other 
fields such as the quantum chromodynamics (QCD) and the quantum hadrondynamics (QHD). Discussions on these 
subjects and practical applications of the coherent-state formalism described in this paper will be anticipated in the 
future. 
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